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ABSTRACT 
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TNiTPnniirTTnKi 



Battery al ignment is one o-f the critical -factors in 
the control o-f the weapons systems on board naval warships. 
The spatial relations o-f the guns, directors, and launchers 
are determined by re-ference to lines and planes -fixed with 
relation to the hull and the determination o-f these 

relations constitutes the battery alignment. This alignment 
only compensates -for the static relations between the 
system elements, and so, when these planes and lines 
distort in a seaway, the alignment is no longer correct. 

Normally, the re-ference lines used -for the alignment 
are the ship centerline -for train and the plane o-f the 
primary element, usually the director, -for elevation. Uhen 
the alignment is done, the ship is alongside a pier in a 
known condition o-f loading and in calm seas. At that time 
the relations, vertical, horizontal, and transverse, 
between the directors and the weapons are measured and 

recorded, and the -fire control system is adjusted so it 
compensates -for static parallax when making the -fire 
control solution. Once the ship goes to sea, these 
relations change as -fuel is burned, stores used and the 
ship moves in waves. The sea, in addition to causing rigid 

body motions o-f the ship hull in heave, pitch, surge, sway, 

roll and yaw, also causes it to -flex in the vertical and 
horizontal planes and to twist about its longitudinal axis. 
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This -flexing distorts the re-ferance lines established in 
the battery alignment and introduces errors in the 
targeting solution arrived at by the -fire control system. 
These errors are due to dynamic parallax, resulting -from 
relative target location due to ship attitude and -flexure. 

Studies by the Naval Sea Systems Command, under 
contract to Rockwell International C13, have shown that the 
-flexure o-f the ship hull o-f monohull vessels appears to be 
s i gn i -f i can t , relative to rigid body motion. It seems that a 
method to measure and compensate -for these errors is 
necessary to improve the -fire control solution and so to 
increase the -fighting e-f -f ec t i veness o-f the ship. 

Measurement o-f the -flexure could be done directly, -for 
example by using a system o-f sighting instruments, and the 
results -fed to the -fire control system as they happen. 
Suitable instruments might be Lasers, light beams or strain 
gauges. All o-f these require the primary equipment, along 
with its support systems, which are subject to -failure and 
in addition the system must be aligned and kept in 
alignment. This results in an i ncreased maintenance load on 
ship's personal, and the reliability problems inherent in 
mechanical devices. In addition, this method does not 
provide an easy capability to predict what the -flexure 
might be at some -future time, thus allowing correction o-f 
the -fire control solution at the actual instant o-f -fire. 

An alternative method involves use o-f a -filter to 
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extract de-flection in-formation -from parameters already 
measured on board the ship. The Shipp's gyroscope provides 
ship heading and speed in-formation and in addition it can 
provide heave, pitch, roll, sway, and yaw in-formation. 
In-formation regarding wind speed and direction is also 
available. Additional in-formation, i -f needed, might be 
•found -from acce 1 erome ters placed at the bow, the stern or 
other critical locations. The advantage o-f using solid 
state devices, instead of directly measuring the 
deflections, is increased reliability and descreased 
maintainace. In addition, the filter formulation lends 
itself to prediction of future distortions, which can be 
fed to the fire control system. 

The objective of this project is the formulation of a 
mathematical model of the ship/sea system that can be used 
with a filter, the formulation of a filter to estimate the 
deflections, and a predictor to estimate future 
distortions. The equations discribing the ship can be 
written as a set of first order differential equations, 
which lend themselves to a state-space model formulation. 
It is then natural to use a Kalman filter for the 
estimator. This method has already been used to estimate 
the rigid-body motions of a ship for input to an automatic 
landing system [2,3,4]. 
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Chk^PTER 1; FORMULATION 



1.1 Literature OvervieiAi 

Wave induced ship uibration is the subject of current 
interest to many investigators. The beginnings of the 
majority of the recent work dates from the late 19<60's, 
when the performance of large bulk carriers was 

investigated by Aertssen [5], Aertssen and De Lembre [6], 

Bell and Taylor C7], and Goodman CSl. These papers were 
mainly concerned with the statistical nature of the 
stresses and bending momients and their prediction, 
Aertssen, and Bell and Taylor doing full scale mesurements 
and Goodman taking a theoretical approach. 

The general trend in the theoretical approaches has 
been to use a modal decomposition, as done for example by 
Bishop and Price, et. al. [9,10,11,12,13], Chen C14], 

Jensen and Pedersen [15], and Kagawa and Fuj i ta [16]. This 

has an advantage when working with regular waves and for 
separating out the contributions of the various modes. 

The problem of the exciting hydrodynamic forces has 

been addressed by Salvesen, Tuck, and Fait insen [17] for 

/ 

the rigid body motions, using strip theory, providing 
linear expressions for the distribution of the wave forces 
over the hull. Skj0rdal and Faltinsen [18] have developed a 
linear theory for the hydrodynami cs of a springing hull. 
Jensen and Pedersen [15] provide a non-linear formulation. 
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Betts, Bishop and Price [1?] tie the strip theory approach 
together with structural representations. 

The use o-f state space techniques -for the prediction 
o-f hull vibration does not seem to have been used by any 
researchers as no re-ferences were -found by the author 
discribing this approach. 

1.2 Ship Mode 1 

There is no established procedure on hoi>j to best model 
the ship and the sea in order to predict -flexural 
distortions. One major consideration is the size o-f the 
models which can make the use o-f state-space techniques 
very cumbersome. For this reason, -for example, it is best 
not to use exact structure and hydrodynamic models, such as 
■finite element techniques, and instead to try to reach a 
compromise between model size and modeling accuracy. A 
question arises when trying to combine a structural model, 
such as modal representation, with a random sea. The sea, 
because o-f dispersion, has a wavelength which is a -function 
o-f wave -frequency. Even though the random sea may be easily 
represented by Fourier techniques at a speci-fic point, its 
representation at some other point in space would be the 
result o-f the propogation o-f an in-finite number of 
frequency components. How to include these phase shift 
effects in a modal analysis is not clear. In addition, the 
computational load may be large due to the number of modes 
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required to accurately represent the quasi-static hull 
distortion, given that the natural -frequencies o-f the hull 
are much higher than the wave -frequencies containing the 
major part o-f the wave energy. 

The approach used here is to approximate the hull as a 
quasistatic non-uni-form beam. Using simple beam theory, one 
starts by considering a short section, dx , o-f the hull <see 
-figure 1.1). The governing equation o-f the beam vibrating 
in a plane o-f symmetry is -found by summing the vertical 
-forces acting at point “a" 

m(x) dx d/dt<dy/dt) = -f<x,t) dx - t)(x,t) * x , t ) -►dV/dx dx 



-f(x,t) dx -► dV/dx dx 



< 1 . 1 > 



or 



m d/dt<dy/dt) = i + V' 



< 1 . 2 ) 



Similary, summing the moments about “a“ one obtains 



M'<x,t) = -V<x,t) 



<1 .3) 



Then by using (1.22) in (1.21) one -finds 



m d/dt(dy/dt) = -f - ' 



(1.4) 




f(Xjt) A X 



F i Qure 1 . 1 
Beam Element 
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From simple, 1 inear beam theory, one has the relation 



M = Ely" 



< 1 .5) 



Using <1.24) in (1.23) one obtains 



m d/dt<dy/dt) = -f - (Ely")" 



< 1 . 6 ) 



which is the governing equation -for vertical vibrations of 
a beam . 

If damping is allowed, then one has 



Once one has the governing equation, the next step is 
to express it in a state-space form. To do this, let xl = y 
and x2 = dy/dt. Then with this coordinate transformation, 
< 1 . 26) becomes 

dxl/dt = x2 



dx2/dt = -<EIxl")" / m- bx2/m+f 

The next step is to represent the space derivative. 
Using a finite difference approach, one lets h = ox, the 



m d/dt<dy/dt) + b dy/dt + <EIy")" = f 



<1 .7) 



< 1 . 8 ) 
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distance between points on the hull. Then i -f there are N 



points, 



L = (N-l)h 



<1.9) 



where L is the length o-f the ship. Expanding the space 
derivative, one has 

<EI xl"")'"' = El'" xl"" + 2 El" xl""" + El xl"""" <1.10) 



Now each part can be represented seperately by -finite 
d i -f -f er ences , with El, El", and El"" as input data. 

At the ends o-f the ship, the conditions o-f zero shear 
and zero bending moment lead to 



El xl"" = 0 
<EI xl"")" = 0 

= El" xl"" -► El xl""" 



< 1 . 11 ) 



Equations <1.29) with boundary conditions <1.291) can 
be represented to the order desired and will result in the 
-formation o-f a set o-f equations, written in matrix -form as 
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dXl/dt = A1»X1 + B1 *E 
Y = C1*X1 



< 1 . 12 ) 



where 



X' = <xl 1 ,xl2,x21 ,x22, . . . ,xNl ,xN2) 
xij = X sub j at point i 

Y = vector o-f responses at the N points. 

A = a 2N X 2N matrix containing the -finite 
di-f-ference representation. 

E' = <-f 1 ,-f2, . . . ,-fN) 

B = distributes the -fj^s to the proper state 



- 14 - 



1.3 Sea Spec tnjm 

Ocean waves are a function of both the wind speed and 
duration, among other parameters. One of the commonly used 
wave spectra that attempt to account for these two factors 
is the Bre tschne i der Spectrum, in which the intensity is 
specified by the significant wave height, H, and the 
duration by the modal wave frequency, Wm. This spectrum has 
the followinq form. 



S<w) = 1 .25 H**2 lAkTi**4 expf-1 . 25< Wm/w) 
4 w**5 



(1.13) 



The movement of a ship through the water causes it to 
experience excitation not at te wave frequency, w, but at 
the encounter frequency, we, where 



we = w - kU cos(beta> 



<1.14) 



The energy of the waves remains the same, however, so 



that 
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S(w)ck«j = S(we>(±Aje 

= S<w) dw dwe 

dwe 



<1.15) 



there-fore 



S<we) = S<w) 

dwe (1.16) 

dw 



where the modulus o-f dwe/dw is used because S<we) is single 
si ded. 

Tr i an ta-fy 1 1 ou and Athans C23 have -found a rational 
approximation to <1.30) in the -following -form 



S<we) 



1 . 25 H*-*2 B<alpha) we**4 

4lAlm < 1 <we/wo) •*•*4) -**3 



<1.17) 



where 
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(1.18) 



alpha = U UAn cos(beta> 
9 

B(alpha) = 1.9339 

1+2 alpha 

wo(alpha) = UJm( 1 + a1 pha) 
0.8409 



U .19) 

U .20) 



The corresponding causal system is 



H(s) = So s**2 

< 1 +2 J( s/wo) + ( s/wo) ** 2 ) **3 



( 1 . 21 ) 



where 



So = 1 . 25 H^-^2 B(alpha) (1.22) 

4 Wm 

J = 0.707 

By writing (1.38) in standard state space -form (see 
chap 2) one has 



dXs/dt = As*Xs + Bs*p 
ys = Cs*Xs 



( 1 .23) 
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I 



where p is white noise o-f intensity q and 



As 



; d 1 8 

! ~wo-**2 “2 Jwo 0 

; 8 0 8 

0 0 -wo**2 

; 0 0 0 

; 0 0 0 



0 


0 


0 


wo*»2 


0 


0 


1 


0 


0 


-2Jwo 


0 


wo**2 


0 


0 


1 


0 


-wo**2 


-2Jwo 



< 



Bs'= [ 0 



0 0 



0 



0 p i ■*wo**33 < 



Cs = [ So 0 



0 



0 



0 0 3 ( 



Since white noise has a constant spectrum tor 
•frequencies 



Sp(w) — q 

TT 



< 



and 



Sp<w) = H<w) **2 S<w) 

= wo**3 So we**4 wo**4 < 

<wo»*4 t we**4)**3 



one may ideniti-fy the intensity as 



q = p i wo-»-*3 So 



< 



• 24a) 

.24b) 

.24c) 

al 1 

.25) 



.26) 



I .27) 
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./■ 



ji 




and the driving noise is taken as 



p ~ N<0, 1 ) 



< 1 .28) 



with pi wo**3 So included in the model B matrix. 
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1.4 Sea Mode 1 



The representation o-f the sea must take account o-f the 
dispersion o-f the waves in space and still allow -for their 
random nature. I -f the relation between driving point 1, pi, 
and distant point 2, p2, is dependent on a -fixed wave 
number, then the sea at p2 will not be properly related to 
pi. To -find a relation that will give the correct relation 
between pi and p2, one starts with Laplace's equation and 
the -free sur-face condition: 

Qxx -^ Qzz = 0 in the -fluid 

< 1 .29) 

Qtt -► g Qz = 0 on the -free sur-face 



Then by assuming that the potential can be written as 



Q<x,z,t) = LKx,t) exp(kz> <1.30) 



and substituting <1.41) into the -free sur-face condition, 
one -finds 

Qtt + gQz = <Utt + gkU) exp<kz> <1.31) 



I -f the deep water dispersion relation, k = w**2/g, is used, 
< 1 . 42) reduces to 
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Utt + w**2 U = 0 



or 

Utt - Utt = 0 



showing that the -free sur-face condition is satisfied as 
expected. Then substituting <1.41) into the Laplace 
equation, one obtains 

g**2 Uxx + Utttt = 0 <1.32) 



It is intresting to note that <1.43) may be obtained 
•from the wave equation by using the expression -for wave 
phase velocity in deep water. 



and 



■from which 
Utt 



c = w / k 



Utt - c**2 Uxx = 0 



w»»2/<w*»2/g) *»2 Uxx 



0 



or 



w»*2 Utt - g**2 Uxx = 0 
and since w*»2 U = -d/dt<dU/dt) 

Utttt g**2 Uxx = 0 



Equation <1.43) can now be represented by -finite 
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di-f-ference equations in space and as a set o-f -first order 
equations in time as 



where 



dXv/dt = Av»Xv + Bv*ys 
E = Cv»Xu Dv*ys 



. <1.33) 



Xv' = <ul 1 ,ul2,ul3,ul4,u21 , . . . ,uN4) 

u j 1 = Uj at the j-th point 

uj2 = dUj/dt 

uj3 = d/dt<dUJ/dt) 

uj4 = d/dt<d/dt<dUj/dt)) 

E = vector o-f wave elevations at the N points 
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1.5 Force Model 



The hydr odynam i c exciting -forces on a ship hull by sea 
waves have been studied for many years. The approach that 
seems the most fruitful for this project is the strip 
theory. Of the several forms available [16,19], the method 
of Oerritsma and Beukelman [19] will be used here. 

The hydrodynamic exciting force on the ship may be 
discribed as being the result of the difference between the 
deflection of the ship hull and the wave profile, z=y-w . 
Then Oerritsma and Beukelman have proposed that the force, 
F<x,t), acting on the ship may be written as 



N(x) * Dz rgB(x)z <1.34) 
Dt 

where 

D = ^ ^ dx = - U 

dt dx dt dt dx 



F( X , t ) = - 



D 

Dt 



I Ma<x) 

L 



Then, using z=y-ui, and writing the -force as 



F<x,t) = -H<x, t) + Z<x, t) 



where 



H<x,t) = Ma<x) D^y 

Dt* 



N<x)-U dMa Dy -*• rgB(x)y 
dx Dt 



Z < X , t ) 



Ma<x) D^w 


N<x)-U dMa 


Dw + rqB<x)w 


Dt^ 


dx 


Dt 



< 1 .35) 



< 1 . 36 ) 
U .37) 



the force can be separated into the part due to the hull 
motions and the part due to the waves. The part due to the 
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hull vibrations may be moved to the le-ft hand side o-f the 



equation and included in the structural model. The wave 
exciting -forces then may be written in a -form that may be 
used in the state-space -formulation. Then the C-matrix and 
D-matrix o-f the sea model are modi-fied with equation <1.53) 
to output the -force instead o-f the elevation. The output 
matrix o-f the spectral model is modi-fied to output the 
-first and second time derivatives o-f the wave elevation in 
addition to the elevation itsel-f. 

1.6 Total mode 1 

The models -for the sea spectrum, the sea sur-face, and 
the ship can be combined into one model by identi-fying the 
output o-f one with the input o-f another in cascade -fashion. 



where 



dX/dt = A*X -t- B*p 
Y = C»X 



< 1 .38) 




Bl*Cv Bl*Dv*Cs 
Av Bv»Cs 

0 As 



B' = C 0 



0 Bs' ] 



C = [ Cs 0 



0 ] 



X^ = < XI ' Xv' 



Xs' ) 



p ~ N<0 , 1 ) 



- 24 - 



CHAPTER 2! STATE SPACE TECHNIQUES 



2.1 State Space Models 

In order to apply Kalman Filtering techniques, one 
must -formulate the problem in the -form o-f a state-space 

model. This is a set o-f -first order d i -f -f er en t i al equations 

that discribe the dynamics o-f the system, its inputs, and 
its outputs. 

The state o-f a system may be de-fined at any time as a 
minimum (not neccessarily unique) set o-f parameters which 
must be speci-fied to di scribe that system. A state vector 

is a vector composed o-f these parameters (state variables). 

The number o-f state variables required is equal to the 
order o-f the system, n. 

For example, a simple mass-spring-damper system is 
discribed by a second order d i -f -f eren t i al equation. 

m*d/dt(dz( t)/dt) + b*dz(t)/dt -► k*z(t) = -f(t) (2.1) 

Then one may de-fine the state variables as the 
position, z, and the velocity, dz/dt, o-f the mass. Next let 
xl = z and x2 = dz/dt. With this change o-f variables, (2.1) 
may be written as 

dxl/dt = x2 



dx2/dt 



-k/m*x 1 



b/m»x2 -► f/m 



( 2 . 2 ) 



or in matrix notation 



where 



dX/dt = A»X + B*F 



A = 

= 

X' = 



: 0 1 ; 

; -k/m -b/m ! 
[ 0 1 ] 

< xl x2 ) 



< 2 . 3 ) 



Al ternat i ve 1 y , because the state variables are not 
unique, one may de-fine the state variables as the -force in 
the spring, K, and the momuntum o-f the mass, P. A 
transformation may be found frc>m one to state vector to the 
other by using 



K = k*z 
P = m*dz/dt 



( 2 . 4 ) 



One may also define an input vector U, containing the 
r inputs to the system and an output vector, Y, containing 
the m outputs. The system equations are then a set of n 
first order differential equations 



dxl/dt = fl<X,U,t) 



< 2 . 5 > 



dxn/dt 



fn<X,U, t) 



or 



dX/dt = F<X,U,t) 


i2.6'f 


The output equations relate the output vector, 
the state vector and the input vector 


Y , to 


yl * gl<X,U,t) 

■ ■ 

a • 

■ ■ 

yiTi = gm < X , U , t ) 

or 

Y = G<X,U,t) 


(2.7) 
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If the system is linear or linearized, one may write 
the functions fi and gi as linear combinations of the 
states and the inputs. 



dx l/dt 




a i J *x j 




b i j *u j 



yk 




ck j «x J 




dk j *u j 



i = [ 1 ,n] 
k = C 1 ,m] 



I n matr i x 



notation this becomes 



( 2 . 8 ) 



(2.9) 



dX/dt = A*X + B*U 
Y = C»X + D*U 



( 2 . 10 ) 



where A i s an n x n matrix di scribing the system 
dynamics, B is n x r and discribes the distribution of 
inputs, C i s an m X n matrix discribing the distribution of 
states to the outputs, and D is a m x r matrix discribing 
the distribution of inputs to outputs. The form discribed 
by (2.10) is known as Standard State-space Form. 



2.2 Kalman Filter 

In this problem the object is to be able to estimate 
the flexural response of the ship hull at time t and from 
this to predict what the response will be at time t+dt . To 
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do this, one common method is to use the Kalman Filter. 
This is a method which derives an optimal (in the least 
squares sense) estimator -for the state o-f a system 
represented by a state-space model . 

Assume one has a mathematical model o-f the system 
under consideration, 

dx/dt = A*x + B*w 

(2.11) 

P = C*x 

where w is white noise o-f intensity Q, the expected value 
o-f X at t=0 , E[x<0)], is x0 = 0, and the covariance o-f x at 
t=0 , EC X < 0 ) X < 0 ) 3 , is X0 . Assume -further that one has 1 
measurements, 2 , which are related to x by 

2 = M*x -►V (2.12) 

where v is measurement noise, modeled as white noise o-f 
intensity R, uncorrelated with the process noise, w. 

Also assume that the system is observable, meaning 
that the initial condition o-f any state may be -found by 
measuring the output -for a -finite period o-f time. 

Then one may propose an observer o-f the -form 

dy/dt = A*y -i- K*< 2 - M*y ) (2.13) 

y(0) = EC x<0) 3 = 0 

K is known as the Kalman gain matrix. 
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4 






I 






*0- 




One may then de-fine the error in the state estimate as 



e = y - X 

and the error covariance as 
P = EC ee' ] 



<2.14) 



<2.15) 



It may be shown C23] that the equation which controls 
propagation o-f the error covariance is 



dP/dt = <A-KM)P P<A-KM)" + KRK' BOB' <2.16) 

P<0) = X<0) =X0 

Examining the covariance matrix, P, -finds that the 
diagonal is made up o-f the variances o-f each o-f the states, 



P = EC ee' 



e 1 
e2 



el e2 ... en ] 



en J 



: el**2 el*e2 . . . : 

: e2*el e2**2 . . . . 

■ • • 

en**2 : 



<2.17) 



So i -f one takes he Trace o-f P 



Tr(P) 



<2.18) 




e i **2 



one has the sum o-f the squares o-f the errors in each o-f the 
states. Athans C231 shows, using the matrix minimum 
prinicipal, that the K that minimizes the Trace o-f P, ie; 
the sum o-f the squares o-f the error, is 



K = PM"R <2.19) 

I -f , as in the present case, none o-f the matrices A,B,M,Q 
and R are -functions o-f time, then dP/dt —^9 in the steady 
state solution, and one has 

0 = AP -»- PA' * BQB' - PM'R MP <2.20) 

Summarizing, then, i -f one has a linear system 
described by 



dx/dt = A*x -► B*w 
p = Cx 



< 2 . 21 ) 



one may estimate, with minimum square error, the state x by 
using y which is described by the equation 



dy/dt = A*y + K<z-My) 



< 2 . 22 ) 




# 



4 



% 



i 







I 







^ I 









with P given by <in the steady state case) 



0 = AP + PA" + BQB" - PM'R MP <2.24) 

Software packages exist, for example CTRLC C26], which will 
solve 2.24 for the steady state values in P and so for the 
values in K . 

Once the values of the state at time t0 are known, 
either from direct measurement or from estimates, the next 
task is to predict the the states at some future time, 
tf=t0+dt. To do this, one may solve 2.291 subject to 
x<t0)=x0 as 



Ncnw, suppose that at time t0, one wants to predict 
x<t0+dt). Then 



but t0<s<=tf is some future time and w<s) is unknown 
because w is random. In fact, since w is white noise and so 
completely random, it is completely unknown and the est 



x(t) = exp<A(t-t0))*x0 + 




(2.25) 



tf 



x<tf) = exp<A<dt) >*x0 + 



exp CA< tf-s) )Bw< s) ds <2.26) 
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estimate o-f it is its average value or expected value, 

EC w<t) ] = 0 <2.27) 

There-fore, the best estimate o-f x(t-f) is 

xp<t-f) = expCAdt>*y(t9) <2.28) 

where y has been used in the initial condition because it 
is obtained -from the Kalman -filter. 
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2.3 Approximations 



The question ot the order o-f the -finite di-f-ference 
representation will be addressed next. The accurarcy o-f the 
distortion results are dependent on that order. For just 
one point o-f the ship model, say j, equation 1.27 can be 
written as (with j understood in the coe-f -f i c i en ts) 



dxjl/dt = xj2 

dxj2/dt = xj\"" - El' xjl'" ~ 

m m 

El " xj!"' -_^xj2 + ^ 

m mm 



(2.29) 



Since the highest derivative in x is the controlling -factor 
in the discretation error, the case o-f a uni-form ship, with 
El' and El" both equal to zero will be investigated. The 
method used will be that used by Smith C23]. 

At point j, (2.30) becomes 



dX/dt = A*X" " + D*X + F 



(2.30) 



where the matrix A is 

0 0 



A = 



and D i s 



D = 



-El/m 0 

0 0 

0 -b/m 
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and F is 



F = 



0 



[ f/m 1 

Now write the tourth derivative o-f X as 



= X<j + 2) -4X<j + l) + 6X<j) - 4X(j-l) + X<j-2) 

h**4 (2.31) 

Using <2.32) in <2.31), one -finds 

dX<j)/dt = A X<J^^2) -4AX<J-U) + 6A X<J) 

h**4 

- 4AX<J-1) + AX<J-2) 
h**4 

D X<j) <2.32) 

and expanding the X<j+1) by Taylor series 

X<j + 1) = X * h*X" 1 h **2*X^^ + 1 h**3*X^"^ -t- ... <2.33) 

2! 3! 

Substituting <2.34) into <2.33), keeping terms to order 
h*»8 and collecting terms leads to 



dX<j)/dt = <^*X" " + D*X F + 

< l*A*X^»h**2 ) + < 1 * A*)<^*h**4 ) 
6 80 



<2.34) 
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where the term in the -first bracket can be interpeted as 
the order h**2 error and that in the second bracket as the 
order h**4 error. 

In order to see the e-f-fect on a particular solution to 
(2.30), consider the solutions o-f the natural modes o-f a 
uni -form beam C22]. The space dependent part is composed o-f 
sinusoids and hyperbolic -functions in the variable kx. 

w(x) = sin(kx) -^ sinh(kx) -^ s*< cos(kx) -^ cosh(kx) ) (2.35) 

where s is a constant dependent on the mode . I -f L is the 
length then 

kL ~ <j-^ )*Pi 

de-fines the approximate values o-f k. Taking the eight 
dirivative o-f (2.35) brings out a k**3 term, so 

w = k**8*w 

and (2.31) can be written as 

dX/dt = D»X F E (2.34) 
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where 



E = _1__ A*X^ *h**4 
80 

= J A*k*»8»X*h»*4 

80 



Then since 



k**4*X = X' 



one may write 



dX/dt = A*X'''' < 1 + J__*h*»4»k*4 ) + D*X + F 

80 

Then the relation between the point spacing and the 
may be found from the second term in the brackets 

e = J «h**4*k*»4 

80 



If the maximum desired error is eo, then 



h**4 < 80 * eo 
k*»4 



<2.37) 



error 
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h < 2.99*L 

< J+i-)Pi 



eo 



(2.38) 



By similar means the relation for order h— squared can 
be found to be 

h < 2 . 45*L eo (2.39) 

( j+i-)Pi 

Table 2.1 shows h vs. eo for both orders for the 
springing mode, and table 2.2 shows h vs. j for a fixed 
eo=0.125. Length is 180 meters. 

Table 2.1 
h vs. eo (j=l) 



h**2 h«*4 



eo 


h 


Np ts 


h 


Np ts 


0.5 


36.76 


3 


53.35 


2 


0.25 


26.00 


4 


44.87 


3 


0.125 


18.36 


6 


37.73 


3 


0.0625 


13.00 


8 


31 .72 


4 



( 



Table 2.2 



hvs. J <eo=0.1 25) 



h **2 



h**4 



j 



h 



1 

2 

3 

4 



18.38 
11 .03 
7.88 
6.13 



Np ts 



h 



6 

9 

13 

17 



37.73 

22.64 

16.17 

12.58 



Np ts 



3 

5 

7 

8 



The most obvious -fact to be gained -from Table 2.1 is 
that -for a given accuracy, the number o-f points and hence 
the number o-f states is much less with the order h**4 
approximation. Each added point requires six additional 
states in the total model , two -for the ship model and -four 
•for the sea model. However, only -five points are required 
to approximate the -fourth direvative to order h»*2 but 
eight points are required -for order h**4. These numbers set 
the minimum number o-f points required -for any accuracy. 

As the number o-f nodes in the vibration increases, the 
number o-f points -for a given accuracy increases rapidly, as 
seen in Table 2.2. Again, the order h**4 approximation is 
the less restrictive o-f the two. It only requires 8 points 
to achieve 12. 5X accuracy with order h»»4, while 17 points 
are required -for the order h**2 approximation. 

This means that while it is possible to use an order 
h**2 approximation, an actual i np 1 emen tat i on will require 
both a higher order and a greater number o-f points. 
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CHAPTER 3: NUMERICAL RESULTS 



3.1 Stability o-f Models 

The mathematical models o-f the sea spectrum and the 
ship structure are stable as modeled. However the force 
model is not. The fourth order equation discribing the sea 
surface has four solutions, two of which are exponentials, 
and it is these solutions that can cause instability. The 
equation is 



Utttt 9**2 Uxx = 0 <3.1) 

which has solutions in the form 
U<x,t) " < exp(iwt), exp<-iwt), exp<wt>, exp<-wt) > 

<3. 2) 

The exponential roots, in the analytical problem, are taken 
care of by the initial conditions. When using a computer, 
however, even if the initial conditions are set equal to 
zero subsequent round-off errors introduce fictitious 
non-zero initial conditions, which excite growing 
exponen t i al s . 

To see this, consider the simple example 

Ytt - a**2 Y = 0 (3.3) 

Y<0) = 0 
Y< infinity) = 0 

which has the analytic solution 

Y(t) = A exp<at> + B exp<-at) (3.4) 
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By appr ox i ma t i n Q this by -finite di-f-ferences, one obtains 



Ytt Y<i-»1) -2*Y(i) Y<i-1) = a**2 Y< i ) <3. 5) 

h**2 

or 

Y<i-H) -<2 <ah)**2> Y< i ) -► Y< i -1 ) = 0 

To solve this -finite di-f-ference equation, take 

Y<n) = r**n 

then one obtains 

r**2 -<2 -► <ah)**2) r -► 1 = 0 
and then -for h < 1/a 

r 1 ah -t- 0<ah**2) 

and 

r ~ 1 - ah -► 0<ah**2) 

Then taking a linear combination o-f the solutions 

Y<n) = A < 1 -*• ah)**n -*• B < 1 ~ ah)*'*n (3.6) 

it can be seen that no matter how small A is, i -f it is not 
exactly equal to zero, the part o-f the solution involving 
the growing exponenitial will eventually dominate the 
solution, swamping the correct solution. I-f the objective 
were only to obtain the solution to the equation, one might 
try integrating backwards in time. Here though, the 

objective is to produce the solution in real time, so the 
equation must be integrated -forward in time. 

This is not as severe a problem as it might seem at -first 
look. One o-f the properties o-f the Kalman Filter is the 
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stab i 1 i zat i on o-f the -filter mode 1 . The ga in that results 
will result in a closed loop observer matrix, 

Acl = A - K»C 

that has all its eigenvalues in the le-ft-hal-f plane. The 
only place the problem arises is in simulation o-f the 
system. By de-fining a reconstruction error, e=x-xhat, and 
using the expressions -for x and xhat, it can be seen that 
the error must obey the di-fterential equation 

de(t) = CA-KC] e<t> (3.7) 
dt 

This means that the r econstruc t i on error goes to zero as t 
becomes large only i -f the closed loop observer matrix is 
asymp tot i cal 1 y stable. 

Simulation o-f the system is necessary to test the 
observer as the model is varied. The method used is to sum 
a number o-f sinusoids, each o-f which has an amplitude 
proportional to the energy in the sea spectra at that wave 
-frequency and a phase that is random with respect to the 
phases o-f the other sinusoids. 

Assume that we are given a sea spectrum, S<w) , and we 
wish to simulate the sea elevation at a number o-f points. 
Since the instantanous sea elevation at a point x and 
-frequency w can be represented as the square root o-f the 
energy at that -frequency (multiplied by 2 since the sea 
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spectra is only defined -from 0 to in-finity), 



V<x,w) = \/2 S<x,w)' 

then the elevation at x, including all -frequencies should 
be 



oo 




<3. 8) 



where p is the random phase. Then including the time -factor 
one has 



CO 




(3.9) 



For a computer simulation, equation (3.9) must be 
discretized to become 



y<x , t) 




S(wj ) ^wj ' 



cost wjt p(wj) > (3.10) 



where here p is a random number uni-formly distributed 
between 0 and two pi . 

The simulated sea sur-face is then used to drive the 
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ship model to -form a simulation model, the output o-f which 
is used as a source o-f the measurements supplied to the 
observer model. This model was run and the results stored 
in a computer -file -for later use. Figure 3.1 shows a 
representative run o-f the simulation. 

3.2 Estimation and Prediction with Correct Model 

The Kalman -filter -for the total model is required to 
construct optimal estimates o-f the states o-f the system 
that can be used to predict the -future output o-f the 
system. This includes those states that cannot be measured 
directly, such as those o-f the sea. To do this it must have 
a source o-f in-formation about the system. 

The -filter gains are a-f-fected not only by the model 
used but also by the noise in the measurements themselves. 
This noise results -from a number o-f sources. For example, 
-from measuring instrument noise, slamming o-f the hull, 
local vibration ef-fects and vibrations excited by external 
sources. The filter then must reject the higher frequency 
noise to obtain the needed information about the states. 
The measurements are found from 

2 = y + V, V ~ N<0,R) (3.11) 

where R is the measurement noise covariance matrix, defined 
as the expected value of the noise process 
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R= E< vv' > 



<3.12) 



Then what is needed are estimates o-f the measurement noises 
that would be encountered in an actual ship hull. Using 
studies o-f ship vibration, the noise covariance . -for 
displacement was taken as 

r<i,i) = 3.0 E -4 meters squared 

at amidships. According to Lewis [27], the amplitude 
amidships will be equal to about 0.25 the value at the end 
■for machinery induced vibration. For velocity o-f 
displacement a covariance o-f 

r<i,i) = 6.0 E -1 meters squared / second squared 

at amidships was used. All o-f-f diagonal values were assumed 
to be equal to zero, implying that the measurement points 
are isolated -from each other. 

The -filter was calculated using the models -for the sea 
spectra, -forces and ship in cascade. This allowed an 
investigation to determine the required measurements to 
insure the observability o-f the system. In addition, the 
■filter pole locations were determined. For the Great Lakes 
Ore Carrier Stewart J. Cort, described in Appendix A, the 
poles were all located at a -frequency o-f less than 1.85 
radians per second. The pole locations are shoiwn in Table 
3.1 -for the output matrix 
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(3.13) 



C 



1 0 0 0 0 0 0 0 
0 1 0 0 0 0 0 0 
0 0 1 0 0 0 0 0 
0 0 0 1 0 0 0 0 



This output matrix provides measurements of the 
displacement. The resulting displacement estimation errors 



are given in Table 3.2 and 


Figures 3.2 
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noisy signal is shotgn in figure 3.4. 
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Figure 3.2 

Estimated vs. Simulated 
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Figure 3.2 
Noisy Measurement 
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Ta^ble 3.1 

Kalman Filter poles 
— ^ 



-6.20E-4 


+ 


i 


2.07 


-7.80E-3 


+ 


i 


1 .08 


-9.80E-4 


+ 


i 


7.38E-1 


-6.07E-3 


+ 


i 


6.46E-1 


-1 .50E-2 


+ 


i 


4.75E-1 


-6.39E-1 


+ 


i 


6.50E-1 


-5.48E-1 


+ 


i 


6.22E-1 


-6.32E-1 


+ 


i 


5.51E-1 


-4.75E-1 


+ 


i 


1 .33E-2 


-7.37E-1 


+ 


i 


4.92E-4 


-6.46E-1 


+ 


i 


4.68E-3 


-5.05E-3 


+ 


i 


4.80E-2 


-6.00E-4 


+ 


i 


4 . SSE-2 



Table 3.2 

RHS Estimation errors 
(meters) 

9.56E-2 
9.36E-2 
1 .54E-1 
9.50E-2 



Prediction o-f the -future values o-f the states, given 
the current measured or estimated values can be done by the 
equation 



dxp = A * xp (3.14) 

dt 

xp(0) = xhat(0) 

I -f A is stable, the predicted value will go to zero as time 
gets large. The objective then is to determine how -far in 
the -future the prediction is acceptable. I -f dt is the 
prediction time step, the actual prediction becomes 
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xp<t) = exp< A *dt > * xhat<t) 



( 3 . 15 ) 



which i nvol ves matr i x-vec tor multiplication only. 

For this system, a number or prediction models were 
tried, but none gave usable predictions. The -first model 
tried was to use the unstabilized system, but as expected, 
the output immediately started growing. 

Next, the system model was stabilized by using a very 
high measurement noise to calculate the -feedback gain. For 
this model, the predicted de -flections showed an initial 
change in the same direction as the true value but rapidly 
decreased in amplitude to about one tenth o-f its amplitude. 
The phase o-f the prediction was very close to that o-f the 
true de-flection. 

Another method tried was to stabilize the model o-f the 
sea sur-face only. This resulted in a stable total model but 
no improvement was seen in the prediction. Here the 
predicted value rapidly diverged -from the true value a-fter 
about -four time steps and the phase shewed no apparent 
relation to that o-f the true phase. 

Finally, a prediction was tried by using only the 
model o-f the ship. This was even less success-ful as the 
prediction immediately diverged -from the true, o-f ten 
increasing when the true value was decreasing and 
vice-versa. 

Though prediction was not acheived by any o-f these 
models, this does not mean that prediction can not be done. 
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It is one o-f the underlying ideas of the state space 
technique that all necessary information about a system is 
included in the state of the system at any time. Prediction 
models that bare investigation are; <1) stabilize the sea 
surface model by treating i t as a regulator problem and 
selectively choosing the elements in the control weighting 
matrix to limit the wave amplitude, <2) try the same method 
on the total model, <3) a useable, if short prediction may 
result from the second method tried if the time step is 
reduced, <4) use the first method tried but scale the 
predicted output. 



- 52 - 



CHAPTER 4: DISCUSSICiN 



The problem addressed in this thesis is the estimation 
of ship hull vibrations -from a minimum number o-f measured 
parameters. These estimations are then used to predict the 
■future de-flections o-f the ship hull. 

The -first step was to develop a mathematical model o-f 
the ship-sea system in state-space -form. This involved 
three separate models. The -first, which models the sea 
power spectrum, was taken -from the literature. This model 
is based on the ITTC ocean spectrum that uses as parameters 
the modal wave -frequency and the signi-ficant wave height. 
In this spectra, the modal -frequency can be taken as 
4.47/Tave, where Tave is the average period measured at the 
calm water level. The signi-ficant wave height may be taken 
from table 4.1 if no measurements are available. 



Tabl e 4.1 

wind speed (knots) 

20 

30 

40 

50 

60 



sig. wave height (m) 
3TI 
5.1 
8.5 
1 1 .0 
14.6 



The model to include the dispersion relation for sea 
waves was developed from a velocity potential 
represen tat i on of the sea, resulting in an unstaola model. 
The instability was due to the increasing exponential 
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solution to the resulting -fourth order in time di -f -ferent i e.1 
equation. 

The third model models the ship hull itsel-f. This was 
derived by assuming that the ship could be represented as a 
simple beam. Then a composite, continuous time, discrete 
space, -finite di-f-ference representation o-f the simple beam 
equation was used -for the -final model . The model was 
developed using a second order -finite di-f-ference 
approximation but -for an actual application, a -fourth order 
approximation should improve the accuracy to a signi-ficant 
degree . 

Gerritsma and Beuke 1 man '' s relation relating the wave 
height to -force on the ship was used to couple the sea 
model with the hull model. Since in a state-space 
-formulation the time derivatives o-f the wave height and the 
de-flection are available in the state vector, implementing 
the -force model simply required additional terms in the 
output matrices o-f the sea spectra and sea sur-face models, 
and added terms in the hull system matrix. 

Using these models, CTRLC C26D, a so-ftware package -for 
control system design, was used on a Digital VAX- 1 1/782 at 
the Joint Computer Facility to do the simulations, develop 
the -filter gains and evaluate the sensitivity to errors in 
the model. The estimation o-f the de-flections o-f the hull 
could be done well with the state space model developed. 
The prediction o-f the de-flection to some time in the -future 
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was 

with 



not successful using the models developed here, 
improved models it could be possible. 



bu t 
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APPENDIX 



PREDICTION RESULTS 

The several methods o-f prediction attempted gave less 
than satis+actory results when compared to the simulated 
signal. Figure A.l shows the prediction that resulted -From 
using the system model stabilized by using a very large 
measurement noise. This resulted in a predicted value that 
was much larger than the simulated value. The same general 
pattern resulted when the measurement noise was varied from 
a low equal to the expected measurement noise to as high as 
could be computed. Nor did changing the starting point ot 
the prediction with relation to the phase o+ the simulated 
signal improve the prediction. 

The next figures show the prediction resulting -from a 
model constructed by stabilizing the sea-sur-face model 
be-fore -forming the total model. In Figure A. 2, the 
predicted signal is very similar to the prediction in 
Figure A.l. I-f just the -first -few time steps are plotted, 
it can be seen that the predicted signal does not track the 
simulated signal even -for a short period o-f time, Figure 
A. 3. 

Finally, Figure A. 4 shows the results o-f using only 
the ship model to do the prediction. As can be seen, this 
prediction has no apparent relation to the simulated 
s i gnal . 
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F i gure A. 1 

Prediction using stabilized total model 
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Figure A . 2 

Prediction using stabilized sea-sur-face model 
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Figure A. 3 

Stabilized sea-sur-f ace , -first 9 time steps o-f prediction 
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F i gure A . 4 

Prediction using only the ship model 
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